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The Narrow Escape problem in a flat cyflndrical microdomain 
with application to diffusion in the synaptic cleft 

Jiirgen Reingruber and David Holcman * 

Abstract 



The mean first passage time (MFPT) for a Brownian particie to reacli a smali target in ceiiuiar 
microdomains is a key parameter for chemical activation. Although asymptotic estimations of the 
MFPT are available for various geometries, these formula cannot be applied to degenerated structures 
' where one dimension of is much smaller compared to the others. Here we study the narrow escape 

, time (NET) problem for a Brownian particle to reach a small target located on the surface of a 

■ flat cylinder, where the cylinder height is comparable to the target size, and much smaller than the 

O ' cylinder radius. When the cylinder is sealed, we estimate the MFPT for a Brownian particle to hit 

a small disk located centrally on the lower surface. For a laterally open cylinder, we estimate the 
conditional probability and the conditional MFPT to reach the small disk before exiting through 
the lateral opening. We apply our results to diflusion in the narrow synaptic cleft, and compute 
the fraction and the mean time for neurotransmitters to find their specific receptors located on the 
postsynaptic terminal. Finally, we confirm our formulas with Brownian simulations. 

G^ ■ 1 Introduction 

The problem of computing the mean first passage time (MFPT) for a Brownian particle to reach a small 
■<!;;j- I target located on a surface of a microdomain, also referred to as the Narrow Escape Time (NET) [36l [14] , 

' is ubiquitous in biophysics and cellular biology because it corresponds to determining the forward binding 

rate of chemical reactions jSHl [21 |3H1 [IH] ■ Applications of the NET ranges from quantitative analysis for 
the resident time of receptors in the postsynaptic density [131 [571 [33 [H] : a fundamental microdomain 
associated to synaptic transmission and plasticity |H!, to scaling laws in physics [7], early steps of viral 
ITJi • infection [TTl [TH [TH] , or the hydrolysis rate of activated phosphodiesterase in rod photoreceptors ^221 [13] • 

I Recent analytical approaches lead to asymptotic formula for the NET in a confined geometry |36l [ST] 

' 1301 129| . For example, in a three dimensional domain of volume V with isoperimetric ratio of order 1, and 

with no bottlenecks, the overall NET to an absorbing circular hole of (dimensionless) radius a centered 
at Xs on the surface is |3U] 

^ ' -^a In a + 0(a In a)) , (1) 



where D is the diffusion constant, L{xs) and N{xs) are the principal curvatures at Xs- In the case of 
a sphere, a precise asymptotic expression with the first three terms was recently obtained in [Bl , where 
the 0(1) term depends on the regular part of the Green's function. The NET computations were further 
generalized to the case of several holes [HI [211 [6], and to stochastic dynamics with a potential well [28ll35]. 

However, the NET formula ([Ij cannot be directly applied to degenerated microdomains where one 
dimension is much smaller than the others. This is for example the case for the synaptic cleft separating 
pre- and post-synaptic neuronal terminals (Fig. [T^) , which can be approximated as a flat cylinder with 
height much smaller compared to its width [3]. Furthermore, in retinal rod photoreceptors sustaining 
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night vision, the outer segment contains thousands of piled flat cyhnders that define the photoresponse 
and the fidehty of the vision under dim hght conditions [331 [T31 [23 ■ 

The goal of this paper is to extend the NET analysis to degenerate domains. More specifically, we 
study the NET of a Brownian particle in a fiat cylinder, where the cylinder height h is much smaller 
compared to the cylinder radius i? (/i <C i?), with a small circular hole of radius a centered on the bottom 
cylinder surface (Fig. [TJd). In the first part, we will analyze the NET to exit the cylindrical domain when 
the boundary is reflecting everywhere except at the small hole, where it is absorbing. Due to the radial 
symmetry, the solution of the the mixed boundary value problem can be expanded in terms of Bessel 
functions. For a flat cylinder with ft, <C i? and a, we find that the NET is given by 

1^1 «o(^) , fR 



^ ^ -yaj ^ 41n _ -3 (2) 

aD ^ 8D \ \a J J ' ^ ' 

where the function oq (^) /\/2 e [0.07, 0.25] is depicted in Fig. [5^. Although we derive (g)) for ft < R, we 
expect that it remains a valid approximation until ft ^ i?, in which case ao(^) ^ j and the leading order 
terms in 1^ and ([T]) coincide. We note that the log-contribution in (HI) comes from the local property of 
the boundary at the hole, whereas in ^ it originates from the degenerated geometry. 

In the second part of the paper, we study a cylinder that is open at lateral boundary, and we present 
asymptotic estimates for the conditional probability p and the conditional mean time Tc that a Brownian 
particle reaches the small hole before leaving the domain through the lateral boundary. For example, 
for a fiat cylinder with h ^ a and R ^ a, when the particle starts at the upper boundary at position 
(r = 0, z = ft) opposite to the small hole, the conditional probability p{0, ft) and the conditional mean 
time rc(0,ft) are (dM]), ^ and (|67l) ) 



T(0,ft) 



1^1 ao{h/a) 
«^ ^o(ff) 



2D 



\/2D ^ ,na. ,„ , , 



,(o,ft) 



R^\n{R/a) "^2h' ' ' r(0, ft) 

1_ r(Oh) 2(ln(W- 

RHniR/a) ^ ' ' 



where a = R/a ^ 1, /3 = ft/a ~ 1 and t(0, ft) is the mean time to reach the small hole when the 
cylinder is closed. These asymptotic expressions can be applied to study diffusion in the synaptic cleft, 
where synaptic transmission depends on neurotransmitters that are released at the presynaptic terminal 
from vesicles and activate receptors located on the opposite post-synaptic neuron (Fig. [1]). The trans- 
mission efficiency depends crucially on the conditional probability for a diffusing neurotransmitter to hit 
a functional receptor before leaving the synaptic cleft. 



2 Mean time to find a small target in a bounded cylindrical 
compartment 

We shall now present our analysis to estimate the MFPT t{x) for a Brownian molecule, initially located 
at position x — (xi, X2, z), to escape a cylinder of radius R and height ft (Fig[T|) through a small circular 
hole of radius a located centrally on the lower surface at z = 0. The cylindrical surface is refiecting, 
except for the small hole where it is absorbing. Due to the radial symmetry, the MFPT is a function of 
the radius r = ^ x\-\- x\ and the height z. Using the small hole radius a, we define the dimensionless 
parameters and variables 

^ ^ R o ^ ^r.^ \^\ o 2 

a a a a a-^ 
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and the scaled MFPT 



which is a solution of 



aD 

w 



■T{r,z) 



D 



T{r,z) 



(4) 



ld_ d_ 

X dx dx 



d_ 

dy 



f(x,y) 
f(x,y) 
f{x,y) 



1 

"M' 

0, y 
0, y 



= 0. 



X 

0, X < 1 

, 1 < X < a 

= 0. 



(5) 



Our goal is to obtain a solution for ([5]) and to clarify its dependency on the parameters a and /3. To 
study the shape of the boundary layer, we note that (O corresponds to a heat equation where the total 
amount of heat produced in is one, independent of a and /3. Furthermore, because the scaled radius 
of the hole through which the heat dissipates is one, it follows that t{x, y) has a finite asymptotic limit 
in the neighborhood of the hole for a — >■ oo and j3 ^ oo. 




(a) 



(b) 



Figure 1: (a) EM picture of a synapse showing a synaptic cleft and the two pre and post-synaptic 
terminals, (b) Schematic representation of a Brownian trajectory in a cylinder of height h and radius R 
with reflecting boundaries, except at the small absorbing disk of radius a <^ R, located centrally on the 
lower surface. We are interested in the case of a flat cylinder with R^ a and h ^ a. 



2.1 Equation for the scaled MFPT f{x,y) 

To derive the solution f(x,y) of (O, we consider the two domains 

Vti — {x\x < 1} and ilo = {x\l < x < a}, 
and obtain for f{x,y) the representation 

{f,;(x,y), (x,2/)efij 
(6) 
To{x,y), {x,y)eno. 

To ensure that T{x^y) is a solution of ([5]) in 57, f{x,y) and the flux ^f(x,?/) have to be continuous at 
X = 1, leading to the conditions 

= To{l,y), 0<y</3 

(7) 

, 0<y</3. 



dfi{x,y) 



dx 



dTo{x,y) 



dx 



Using a separation of variable method, we expand Ti(x, y) and ro(x, y) in series 

oo 

n{x,y) = bnqn{x)vn{y) + Wi{x, y) 

ri=0 
oo 

{x)un{y) +Wo(x,y) 

where 



(8) 
(9) 



k -— I 



(n+ i)7r 



Mo = , M„(y) = cos (fc„y) (rt > 1) , t)„(2/) = sin(;„j/) (rt > 0) , 



(10) 



Wi{x,y) and Wo{x,y) are the inhomogeneous solutions of ([5]) that vanish at x = 1, and Pn(a:) and (7n(2^) 
will be derived below in terms of the modified Bessel functions h){x) and _ftro(.x), and are normalized such 
that pn{l) = q„(l) = 1. 

The functions Vn{y) and w„(j/) satisfy the orthogonality relations 



Uniy)Um{y)dy ^ I Vn{y)v,n{y)dy ^ ^5nm 



Vn{y)Um{y)dy = -inm , 



where 



2 



m > 1 



V2 V2 1 



/3Z„ TT 71 + 



— , m = 



is an orthogonal matrix satisfying 



Using the orthogonality relations, we obtain the expansions 

oo oo 



(11) 



(12) 



(13) 



2.1.1 Derivation of p„(x) and Wo{x,y) 
The equation for To{x,y) in i7o is 



1^5 5' ^ " / ^ 1 / ^ o 

9 d 
—fo{x,y) =0, —fo{x,y) =0 

ay y=0,l3 ox x=a 



and we choose Wo{x, y) to satisfy 



I d d 92 



a: 9a; dx dy'^ 
d 

Wo(l,2/)=0, —Wo{x,y) 



Wo(x,y) 



1 



{x,y) £ no 



= 0, —Wo{x,y) 

y=0,p ox 



0. 



(14) 
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The solution for Wo(x,y) is 



lux X ~ 1 



27r/3 4|r2| 

Inserting To{x,y) from ^ into yields for Pn{x) the equations 

-— Xt; fc„ p„(x)=:0, l<x<a 

X ax ox ' 



(15) 



(9 



0. 



Using the modified Bessel functions Io{x) and -f4ro(x) and the relations [5] (/o(x) = /i(x), Kq{x) = 
—Ki{x), we obtain 



Pn(x) 



(16) 



with 



Fo{x,y) - /o(x)A'i(2/) +i^o(a^)/l(2/)■ 



2.1.2 Derivation of qn{x) and Wi{x,y) 

Proceeding similarly to the previous paragraph, the equation for fi{x,y) in fli is 



^ d d 52 \ ^ ^ 1 , , ^ 
d 

f^(x, 0) = 0, —fi{x,y) 



(17) 



0. 



We choose ?«i(x, y) to satisfy 



15 9 92 
X dx dx dy^ 



i{x,y) = - 
d' 



Wj(l,y)=0, Wj(x, 0) = 0, —Wi{x,y) 

oy 



{x,y) £ 

= 0, 



V=I3 



with solution 



00 

n=l 



cosh(z„(/3 - j/)) _ x^ - 1 
cosh(z„^) 
cosh(2;„(/3 - y)) 



(18) 



cosh(z„/3) 

where z„ are the positive zeros of the Bessel function Jo(x), and the coefficients c„ are given by 



JoiZn)'^ Jo 



Jo(z„x) 



x^-1 



-xdx . 



(19) 



To derive expression we used the orthogonality relation [S] 

Jo{ZnX)jQ{ZjnX)xdx = J„™ ^ (-'"o(^n)^ + Joi^nf) 

Inserting Ti{x,y) from ([5]) into (jl7p gives for qn{x) the equation 

'i_a_ _9_ 

X 9x 9x 



x-^ - ll ] q„{x) = 0, X < 1 



'Zn(l) = 1, 



and the solution that is regular at x = is 



qn{x) 



(20) 
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2.1.3 General expression for f{x,y) 

Using the expressions for pn{x), qn{x), Wi{x^y) and Wo{x,y), the NET sohition is 



f{x,y) = < 



To{x,y) = ^ a„ 



Fo(A;„x, fc„a) \nx x^-\ 



-Fb(fcn, fcna) 



27r/3 4|17| 



where the unknown coefficients a„ and 5„ wiU be determined by patching the two expressions at 
The continuity condition for T{x,y) at x = 1 gives 

CO 00 

anUniy) = ^ bnVniy) , 
n=0 n=0 

and using the expansions in (jl3p . we obtain that a„ and 6„ are related by 

00 00 



m— 



n=0 



The continuity condition for the flux at x = 1 gives 
h{ln) / X , FiiK^Ka) 



S^, , h{ln) r \ 1 PiV^n.Ka) , 1 1 



with 



This can be rewritten as 



KPnVn{y) + ^ a„Q!„W„(y) = ^ -fnUniy) ■ 



n=0 



n=0 



n=0 



where 



' Fo{kn,kna) 

and the 7„ are implicitly defined by the equation 



By using the expansions in (1131) we obtain from 

oo oo 

m—0 771=0 

oo 

Finally, using the relations between and bn given in (|22p. we obtain the matrix equations 

oo oo 

m—0 m—0 

oo 



For given a and /3, by truncating and numerically solving these equations we find approximated values 
for a„ and 6„, and from this we obtain an approximation for T{x,y). We will analyze the equations for 
a„ and 6„ in more detail later on. 



2.2 MFPT with a uniform initial distribution 



We shall first consider the average MFPT t{x) when the Brownian particle is initially uniformly dis- 
tributed at radial position x. Using (|2ip we obtain 



f{x) 



T{x,y)dy 



To{x) 



^ oo 



1 



fif^,ln h{ln) 13 Jo 
2-1 



ao In a; 



w^{x,y)dy, x & Qi 



X G rto- 



(28) 



Expression shows that is the averaged MFPT for Brownian particles that are initially uniformly 
distributed at x = 1. The expression for fo(x) has an intuitive interpretation; the escape time starting 
at a; > 1 is the sum of the average time to reach x = 1, plus the escape time starting at x = 1. 
The average MFPT f for particles that are initially uniformly distributed in f2 is 



1 

\n\ 



t{x, y)dV - 



1 



\n\ \n„ 



T{x,y)dV 



1 



(29) 



where 







Wi{x, y)xdxdy 







V2 



4 1 

41na — 3H tt j 



(30) 
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The time fi is the average MFPT for particles starting uniformly distributed in the inner cylinder fli, 
and To is the average MFPT for particles starting uniformly distributed in the annulus flo- We shall now 
derive asymptotic limits for fi and Tq under various conditions. 



2.3 Asymptotic expressions for a cylinder with R ^ a, and a flat cylinder 
with i? > a and h <t: R 

We will first derive asymptotic expressions for f for a cylinder with i? ^ a (a ^ 1) and arbitrary height 
h, and we will then focus on a flat cylinder with h ^ R {/3 ^ a). We show that ao{a,f3) is the leading 
order contribution to f for a 3> 1. For a flat cylinder with a 3> 1 and 13/ a ^ 1, we further have that 
ao{a, f3) « ao{f3). To derive ao(/3) as a function of /?, we consider the limit a — ?> oo while /? stays bounded 
(i? — > oo with flnite h). We show that t{x, y) and f have flnite asymptotic limits for a — > oo that depend 
only on /3. 

We start by considering the limit a ^ 1. The function Wi{x,y) in (fT8)) is of the order ~ 
and can be neglected, and we have 

nix,y) ~ y!^» '^°''.^"^? «n(^). (31) 

Because the average time starting uniformly distributed in fli is similar to the the average time 
f(l) — starting uniformly distributed at a; = 1, the contribution of in (I^Hl) is by a factor 
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smaller compared to the contribution of fo, and we arrive at the asymptotic expression 

00(0:, /?) 4 In a — 3 



The dimensional time r is 



V2 



\V\^ |V|ao(a,/3) , 



aD aD ^ 8D 
In particular, for /3 ^ Ina and a ^ 1 wc obtain the result 



87r/3 



4 In 



- 3 



V2 



\V\ aoia,P) 
aD V2 



(32) 



(33) 



(34) 



Equations show that ao(a, f3) is the leading term that determines the average MFPT for a >• 1. 

To further evaluate r, we shall now estimate ao{a, (3) for a fiat cylinder with a small hole, when /3 <C a 
and a ^ 1, by considering the limit a — 00 while /3 remains finite (i? — )■ 00 with fixed h). For a — 00, 
the scaled times f (a:, y) and f have have finite limits that depend on /3, and only the dimensional times 
T(r, z) and r diverge ~ B? . In this limit, the coefficients a„, /3„ and 7„ in (p4|) and (|25|) are given by 



ao = 0, a 



n — "'n 



'i^o(fcn) 

with 7o = , and (|26l) simplifies to 



j^i(fc„) ^ . , /i(/„) 



ML) 



SnO 



(35) 



5Z (/3n + ar„)e Cn070 ■ 



(36) 



a„, /3„ and 7„ are functions of /3 only, and hence, also a„ and &„ depend only on /3. For a 3> 1 and 
/3/a ^ 1, f(a;, y) in (pi]) is in first order given by 



-(x,2/) < 



ao(/3) 



I V2 '5/"(^)i.o(M 



Koiknx) , , In a; 

Wn(2/) + ^, l<a;<a, 



where we used 



FQ{k„x,k„a) Ko{knx) 



a 3> 1 and x a . 



Fo{kn,kna) Ka{kn) 
We conclude that the NET for a ^ 1 and /3/q! ^ 1 is in leading order 



\V\aoiP) , R 



+ 



4 In 



aD V2 81) 

In the next section we shall analyze the behavior of oq as a function of j3. 



(37) 



(38) 



2.3.1 Behavior of ao{/3) as a function of /3 

To evaluate ao(/3) as a function of /3, we solve numerically by truncating the series at sufficiently high 
values n: in Fig. [2l^a) we plot the analytic approach result for ao{(3)/^/2 and confirm that it agrees well 
with results from Brownian simulations that were performed with a — 50. Interestingly, Fig. [2^ shows 
that the simulation result for (3 = AO (when 13 is comparable to a = 50) still agrees very well with the 
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analytic result derived with the assumption a 3> /?, suggesting that ao(/3) remains a good approximation 
until values (3 ^ a {h ^ R). As a consequence, this suggests that (1551) is an acceptable approximation 
for r until values h ^ R. Fig. shows that approaches the value j for large /3, thus, from 

we recover the narrow escape formula r w |36[ I10[ 131) derived for a volume with isoperimetric ratio 

of order 1. Conversely, ([5^ shows that the validity of the narrow escape formula is not limited to 
the range where h ^ R, but it is already a valid approximation when \na /3 and (3 > 40 (Fig. [2ji). 
Hence, we conclude that r = is a good approximation even for an oblate volume with R ^ h (and 
h :$> a). In the opposite limit /3 ^ 0, we find that ao(/?)/\/2 docs not converge towards zero (Fig. [5^), but 
lim^_5.o ao(/3)/-\/2 « 0.071 (in appendix [K\ we derive an analytical approximation for ao(/3) for /3 — >■ 0). In 
Fig. [2b we show the effect of the truncation level n on the value of ao(/3)/-\/2: after n ^ 100 the steady 
state regime is achieved. Finally, in Fig.[2};-d, we compare the values of the coefficients ao(/3) and bo{(3), 
where 6o(/?) is obtained using ((22|). The graphs show that bo{/3) « ao{(3) is a good approximation, and 
we will use this in section [2.3.31 




Figure 2: (a) Graph of """^ as a function of /3 for a ^ 1. The values for qq are obtained by truncating 

and numerically solving ([36l) (for /3 > 10 we truncate at n = 500). ao(/3)/\/2 is the average MFPT for 
Brownian particles starting uniformly distributed at a; = 1. The analytic estimations for ao(/3)/\/2 are 
compared to results from Brownian simulations (data points) obtained for 10 000 Brownian trajectories, 
starting uniformly distributed at a; = 1 ( for a = 50). (b) Dependency of on the truncation level 

n for various (5. (c) Comparison between ao(/3) and &o(/3), where bo is obtained from ([22[) . showing that 
ao(/3) « boiP). (d) Relative difference (ao(/3) - 6o(/3))/ao(/3)- 
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2.3.2 Boundary layer analysis: particles starting near the absorbing hole 

In the neighborhood of the small absorbing window (for x ^ I and y ^ 1), there is a boundary layer (BL) 
where the behavior of t{x, y) is very different compared to large x and y. In Fig. [3] we study numerically 
the shape of the BL using ((37)) . The different panels depict T{x,y) in the neighborhood of the absorbing 
window for various /3. The plots show that a boundary layer starts to evolve around /3 ~ 0.3, and the 
evolution is almost finished for /3 ^ 10 (there is no significant difference between the plots for f3 — 10 and 
/3 = 50). Furthermore, the approximate extent of the boundary layer for large /S is Ax ~ Ay ^ 10. 

In Fig. 3^, we show t{x, /3) for particles released on the upper surface at y — (3 as a, function of 
X and for various (3. Such a situation is relevant at synapses where neurotransmitters are released at 
the presynaptic terminal, located opposite to the surface with the postsynaptic density (PSD) where 
receptors are clustered [20l |4]. The hole radius a would correspond to the radius of the PSD. Fig. IHa) 
shows that when the height of the synaptic cleft is comparable to the radius of the PSD (/3 ^ 1), t{x,/3) 
changes considerably as a function of the radial release position x. In contrast, for /3 ^ 1 the release site 
is outside the boundary layer and the NET is almost independent of x and is well approximated by 1/4. 




2.3.3 Impact of truncating the series for t{x, y) in (j37p 

We now study the error induced by truncating the sum in p7|) at levels n ~ 1 by considering the truncated 



series 



< 1 



«o(/3) , ,a\^o{hx) \nx 



(39) 



I V2 



1=1 



Koih) 



2tt(3 ' 



1 < X < a, 
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(a) (b) 



Figure 4: (a) The NET t{x, fi) computed from ([37)1 for particles that are released at the upper cylinder 
surface at radial position x for various /?. We further compare the analytic computations with results 
from Brownian simulations obtained with 10000 particles and a = 50. (b) Comparison of f(0,/3) with 
the approximation f(°)(0, /3)ap ~ aQ{(3) / Io{jg) given in PO]) . 



To evaluate the error induced by the truncation, we first compute the coefficients a„ and 6„ with high 
precision (using a truncation level n ~ 200), and then use these values in In Fig. [5l we show the 

effect of the truncations for various n and /?: interestingly, the numerical analysis reveals that for /3 ^ 1, 
truncating at n = or n = 1 already provides a very good approximation. The accuracy of the truncation 
depends on j3, and n has to be increased for larger /? in order to maintain a similar accuracy (Fig. (St). 
In Fig. H^-c, we plot the effect of the truncation as a function oi x ior y = j3 (particles are released at the 
upper surface), and in Fig. |3Ji-f, the starting position y is reduced to y = 0.7/3, y = 0.4/3 and y = 0.1/3. 

Due to the truncation, at the patching boundary x ~ 1, (p9)) has a small discontinuity A'"^(y) — 
f (")(!+ _ f(")(l^,y). For example, for n = and y = (3 (see Fig. [S^-c) we obtain A'^''^(/3) = 
^ - 6o(/3) « ao{l3){l - ^), where we used ao(/3) « &o(/3) (see Fig.[l:-d). 

Finally, for /3 < 1, when a Brownian particle is released at the center of the upper surface {x — O^y — 
/3), using the truncation n = 0, we obtain from (I39p 



(0, /3) « f (0) (0, /?) = ^ - fj^ = (1) , (40) 

^o(.2^) -'ol2:a^ ^o{2p) 



where we additionally used that ao(/3 « 6o(/3) and f(l) = ao/V2. In Fig. |4}d, we test this approximation 
as a function of /3 by comparing it with f(0,/3), computed with high accuracy {n ^ 200). We find that 
this approximation is valid until /3 ~ 1. 

2.3.4 Analogy with the electrified disk problem 

For /3 ^ 1 and a 3> 1, an asymptotic solution for t{x, y) can be obtained by considering the analogy with 
the electrified disk problem: the total outflux —1 through the hole leads to the electrified disk problem 
with a disk charge Q = — 1, and using the capacitance C = 4 of the unit disk [33l [32], we find that 
the disk potential is Q/C = —1/4. Using the solution U{x,y) of the electrified disk problem with disk 
potential —1/4 jSSj . we obtain the asymptotic correspondence 

T{x,y) = U{x,y) + ^, a>l,/3>l. (41) 



Hence, for large /3, U{x,y) determines the shape of the boundary layer. Furthermore, by comparing (j4T 
with (I37p . we recover that ^ = i for ;3 ~ a > 1. 
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(d) 13 = 1 and y = 0.7 X 13 (e) /3 = 1 and j/ = 0.4 X /3 (f) /3 = 1 and y = 0.1 X /3 

Figure 5: The truncated NET f'-"'\x,y) fr om ((39|) for various truncation levels n and several /3. The 
figure shows that for /3 1 very good results can be obtained by considering a low truncation level n 2. 
(a)-(c): T^"\x,P) for various (3. (d)-(f): f^"\x,y) for /3 = 1 and various y. 



3 Conditional probability to reach the small target before leav- 
ing a laterally open cylinder 

When the cylinder is open at the lateral boundary, we shall now compute the conditional probability 
p{x,y) that a Brownian particle, initially at position {x,y), reaches the small target disk before leaving 
the cylinder through the lateral opening. Because the geometry of the synaptic cleft can be approximated 
by a laterally open cylinder [TJ [22l |34] , we will use our computations to estimate the efficiency of receptor 
activation at a synapse. Indeed, at the presynaptic site, vesicles release neurotransmitters into the 
synaptic cleft, and the diffusing neurotransmitter either bind to and thereby activate receptors located 
on postsynaptic terminal, or they leave the synaptic cleft without activating a receptor. We will first 
derive a general expression for the conditional probability p{x, y) to hit the small target before exiting, 
and then compute average values for uniform initial distributions. Finally, we will determine the leading 
order behavior in a flat cylinder with a and h ^ R. 

The conditional probability p{x, y) satisfies the Laplace equation [TB] 

^ d d . . ^ 

^x ox ox oy ) 

'p(a;,0) = 1 forx< 1, p(2;,a)=0 (42) 

d d 
—p{x, 0) = for 1 < X < a , —p{x, f3) — . 
oy Oy 

Similarly to the analysis of T{x,y) in the previous section, we solve ([S]) in the subdomains fli and fio, 
and then patch the two solutions Pi{x,y) and Po{x,y) at the boundary x — 1. The general expressions 
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are 



p{x,y) = < 



n=0 



(43) 



In ^ 



^/ ^"^^ ^ " Go{kn,kna) 



Uniy), l<x<a 



where 



Go{x,y) = Io{x)Ko{y) - Ko{x)Io{y), 



and the unknown coefficients and 6^^ are functions of a and /3 and are related by a relation similar to 
The coefficient resp. 6^ are given by 



m=0 

CO 

^ ^ (/3m + Q!^i)CTrm^n 



C»o7o 



(44) 



m=0 



with 



and 



In a 



Go[kn,kna) hyln) Ina 



(45) 



Gi{x,y) = — Go(a;,y) =/i(a;)ifo(y)+J^i(a;)/o(y). 
For the /?„ we omitted the superscript ^ because they coincide with the /?„ already defined in ([24| . 

3.1 Conditional probabilities with uniform initial distributions 

The fraction p{x) of Brownian particles that eventually reach the target starting initially uniformly 
distributed at x is 



1 

p{x) = ^ J Pi.x,y)dy = < 



1-E 



hi IpjlnX) 
pin UL) ' 



X < 1 



(46) 



V V 2 / in a In a 

where p(l) — 1 — Oq/v^- For particles starting initially uniformly distributed in f2, the average probability 



IS 



p = p(x)xdx — p[l) — 
a Jo 



21na-l 1 



1 2 \ ^ 6n /l(?„) 

2a2 In a ^ ' ^^1^1^ 



(47) 
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When the particles are initially uniformly released at the upper surface within an area x < xq, the fraction 
Pisixo) that will reach the target before escaping through the lateral opening is 



p{x, /3)xdx 



h{ln) 



Xo < 1 



V2J 41na 



2 2 

00 



(48) 



_2_ v^, ,x„<c< / XQGi{knXo,kna) 
4^} kl V Gi(fc„,fc„a) 



1 



1 < xq < a . 



3.2 Asymptotic expressions for the conditional probability in a flat cylinder 
with a ^ 1 and /3 ^ a 

To obtain an asymptotic expression for the conditional probability p in the limit a ^ 1 and (3 ^ a, we 
estimate af^ and b^. For a » 1 and f3/a <^ 1 we have 



In, 



' Ko{k„ 



an, [n > 1) , 



where a„ are given in (I55|) . For the scaled coefficients aj^ = ^^a^ we obtain from the asymptotic 
equations 



1 °° 1 



(49) 



When a is large such that ^ /?„, reduces to ([55t for the coefficients a„ of f{x, y), hence, w a„. 
Because /?„ are monotonically increasing with n, if ^ /3„ for n = 0, this is also valid for n > 0. Hence, 

by setting n = 0, we obtain the condition ^ ^ /o(^) ' ^^^isfied for small (3 when Ina ^ /3, and for 

large (3 when Ina^ (3^. In this limit, the asymptotic expressions for aj^ and b^^ are 



a„ and Oj; = 6„ 

ma ma 



(50) 



Finally, by inserting and 6fj into P5|) we obtain for for In a ^ ^ ^ j the asymptotic formula 

277/3 ^ /o(^«a^) . ^ , , 

Ina'^ /o((„) 



p{x,y) = < 



n=0 



.nMll 2^/3^ j^o(fcnX) 
ma ma ivo(Kn) 



(51) 



1 < X < a 



where p(l) = 1 - 1 - |^|-^ (glD and we used 

GgjknX, kna) ^ KpjknX) 

Go(fc„,fc„a) Ko{kn) 
By considering p7p for f (a;, y), this can be rewritten as 

2ttP 



a ^ 1 and x a . 



v{x,y) 



l-—^f{x,y), x<l 
m a 



m a \ in a / ma 



(52) 
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showing an interesting relation between the conditional probability and the MFPT. As an example, for 
particles that start uniformly distributed at a; = 1 we have 



In a 



(53) 



and using that f(l) ^, we find that the probability p[l) approaches one for large a. In contrast, for 
particles that start initially uniformly distributed in Q, we obtain from (|47p 



2 ma 



(54) 



which shows that p tends to zero for large a, in contrast to p(l). 
3.3 Numerical evaluations 

In Fig. we plot p{\) = 1 ^ as a function of a for various /3 (we numerically solve (|44l) with high 
accuracy up to n ~ 400). As a — >■ 1, the probability p(l) tends to zero because the particles start close 
to the lateral boundary. The asymptotic limit of p{l) for large a is one, as shown in ((53|). However, the 
convergence is only logarithmical and therefore very slow. In Fig. [61d we display p{l) as a function of /? 
for various a: for a fixed value a > 1, the asymptotic limit of p(l) for /3 — > is one, and zero for /3 — > oo. 
The limit for /3 — is intuitive, because the particles start next to the target. To obtain the asymptotic 
limit for /? — )■ cx) we consider (|44l) for aj^ and (|45|) : because the coefficients (n > 1) and (n > 0) 
tend to zero for large /3, the non- vanishing part of (|44| is aQ^noao = 7o?n,o, from which it follows that 
the asymptotic limit of Oq for large /3 is v^, and thus p{l) = 1 — Oq/v^ converges to zero in this limit. 
Indeed, for fixed a and increasing (3, it becomes less probable that a particle reaches the disk before the 
lateral boundary. 
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Figure 6: Average conditional probability p(l) = 1 — for Brownian particles starting uniformly 
distributed at a; = 1 to reach the target before leaving the cylinder through the lateral boundary. The 
coefficient Qq are obtained by truncating and numerically solving (1441) with high accuracy, (a) p{l) as a 
function of a for different /3. (b) p{l) as a function of /3 for different a. 



We now study the conditional probability for particles that are released at the upper surface, which 
is relevant for synaptic transmission. In Fig. [T^-c we plot p{x,l3) as a function of x for a = 10, a = 5 
and a = 100, and various (3 between 0.1 and 10, and in Fig. [TK we further show that our analytical 
computations agree with results from Brownian simulations. In Fig. [7Jl we further depict the average 
probability pjs^xo) when particles are released at the upper surface within an area of radius xq for similar 
values a and /3 as in panel (c). In general. Fig. [7] shows that p{x,j3) is very sensitive to the cylinder 
height P and the release position x. For very fiat cylinders with /? ^ 1, when the particles are released in 
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Figure 7: Panels a-c: Conditional probability p(x,/3) (from (|43)) ) for a Brownian particle to reach the 
target at {x < l,y — 0) when released at the upper surface a.t y — (3 and radial position x. The data 
points in (a) are obtained from Brownian simulations with 10^ particles, (c) The values for /3 are as 
in panel (a), (d) Average conditional probability ppixo) (from (05])) when the particles are released 
uniformly on the upper surface in an area of radius xq with values /? as in (c). 



the area opposite to the hole (for x < 1), the probability is ^ 1, whereas when they are released outside 
this region (for a; > 1), the probability that they reach the target before exiting decreases considerably 
(Fig.[7^-c). For example, Fig.[7)D, obtained for /3 = 0.5 and a = 5, shows that the conditional probability 
is larger than 90% when the particles are released at the upper surface within the area x < 1, whereas it 
decreases to around 60% when they are released at a; ^ 2. For cylinder with large a the impact of the 
release position is much less pronounced (panel (d) with a = 100). In addition, Fig. [7] shows that the 
conditional probability to reach the target is more sensitive to changes in the cylinder height /3 compared 
to the radius a. 

3.4 Impact of truncating the series for p{x, y) in ( 1431) 

We now proceed similarly to section [2. 3.31 and estimate the error induced by truncating the series in (|43l) 
at small n. We first compute the coefficients af and 6f using (|44|) with a high accuracy, and then use 
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these coefRcients to define the truncated probabihty 



l.±tf-f^v.iy), .<1 



1 - 



^/2/ hi 



(55) 



In Fig.[Sl we plot /3) for various a and /3 as a function of x, and show that truncating at n = 0, 1, 2 

aheady gives good approximations (we plot p'"^ (x, (3) for n = 100 to show the error induced by the low 
truncations). Similarly as in section [5331 conclude that truncation at n = 1 or n = 2 already provides 
a good approximation for /3 ^ 1. 




4 Conditional mean time to reach a small target before escaping 
a laterally open cylinder 

After having estimated the conditional probability p{x,y), we are now in a position to study the con- 
ditional MFPT fc{x,y) to reach a target before escaping through the lateral cylinder boundary. This 
analysis will estimate the time scale of the synaptic response. For example, the conditional time Tc{0,/3) 
in (I67p provides an estimate for the mean time until postsynaptic receptors become activated after trans- 
mitter release into the synaptic cleft at the upper cylinder surface. In a wide cylinder with a 3> 1, (|67p 
shows that the conditional time fc(0, f3) is roughly by a factor l/(lnQ:)^ faster than the mean time f(0, /3) 
for a laterally closed cylinder, because in a closed cylinder all trajectories that are reflected at the lateral 
boundary contribute to the mean time, thereby increasing the mean compared to an open cylinder. Thus, 
to obtain a realistic estimation of the synaptic activation time, we have to consider the conditional time 
Tc for an open cylinder. 

To determine the mean conditional time Tc(x, y), we use the known conditional probability p{x, y) to 
solve an equation for the function A{x,y) — Tc{x,y)p{x,y), and then we obtain Tc{x,y) from Tc{x,y) = 
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The function A{x,y) satisfies 0123 HI] 



15 9 92 \^ ^ 
X ox ox oy^ J 

A{x,y) 

dy 



p{x,y) 



ix,y) € n 



, y = ,x < I 
, X = a 



(56) 



, y = 0, X > I 

0, 

Proceeding as in the previous sections, we expand A{x, y) as 



f:b^^-f^v^iy) + wfix,y), x<l 



A(a;,y) = < 



(57) 



^2 Ina ' ^ " Go{kn,kna) 



E 



,(y)+u;^(a:,y), 1 < a; < 



where y) and w^{x, y) are the inhomogeneous solutions of (j56p that vanish at a; = 1. The coefficients 



a„ and b„ are related and 6„ in ([22j), and a„ satisfies 



A 

ni ? 



m=0 



m— 



where 



Ina Go(K„,fe„a) ^o(/„) 



and the 7„ are implicitly defined through 



d_ 

dx 



{w^{x,y) - wf{x,y)) =S^^^Un{y) ■ 



(58) 



(59) 



(60) 



To determine the coefficients 7;^, we first evaluate w^(a;, y) and w^{x, y). When a 1, wf{x, y) is of the 
order a"^ and we neglect its contribution in first approximation. Using p(x,y) from ()43p . the equation 
for w^{x, y) is 



19 9 92 
X dx dx dy^ 



d_ 

dy 



Wo{x,y) 



lii \^ Ina ^ Go{kn,kna) 



, X = 1 and x = a 



(y) 



(61) 



w^(a;, y) = 0, y = /? and ?/ = . 



To solve (|6T|) we expand w-^{x, y) in terms of Wn(y), 

00 

u;„^(x,2/) = ^„^(°)(a:) + J] (x)u„(2/) , 

n=l 

and inserting this expansion into (j6ip gives for w^'"'^^ {x) the solution 



+ 1) - (a^ - Ina - 1)|^ - (Ina + 1) 
ma 



(62) 
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The higher order functions w'^^"\x) (n > 1) satisfy the equation 



1 < X < a 



To proceed, we now truncate the series for 'w^{x,y) in (j62p at n = and use only the first order 
approximation w^{x,y) ~ w^i^^°'(a;). We expect that this aheady provides a good approximation because 
the coefficients are small for for n > 1 and large a, and, as shown in section [3.41 truncation at n = 
already gives a very good approximation for p{x^ y when /3 ^ 1. Hence, we expect that our analysis is a 
valid approximation for large a and small (3^1. With truncation at n = 0, the parameters 7^ defined 
in (l60l) are 



7,f=7^^nO, with 7^= Ml) 



- 2(lnQ!)2 - 21na - 1 
nal^ 



1 



y27r/3 2(lna)2 



(63) 



2^ /q(_il^ 

Similar to (|50p , for In a 3> 7~FT ' ^^^^ asymptotic solutions for and 6^ in terms of a„ and 6„ are 

TT /l(2^) 



" 2(lna)2 
Finally, we obtain for Tc{x,y) the approximation 



and = —Tr^-Krrbr, 



2{\na) 



(64) 



fcix, y) 



p{x,y) p{x,y) 



^,aIo(Ux) p{l) 
2^^^nj^v„{y)^,^^j^Tix,y), x<l 



^2 Ina 



E 



2(lna)2 
j^KoiKx) 



(65) 



ifo(fcn) 



where p{x,y) is given in (|5ip . To estimate how much -fc(a;, y) is faster compared to T{x,y), we consider 
particles that are released centrally on the upper surface at position (0,/3). By taking into account (|52|) 
for p{x, y), we obtain 



fc(0,/3) 



(66) 



Furthermore, using from (|40l) the approximation f(0,/3) 
approximation 



- -^f(l), we obtain for ^ < 1 the 



2/3 ^ 



rc(0,/3) 



1 - 



^/o(^)f(0,/3) f(o,;3) 
l-f^f(0,/3) 2(lna)2 



(67) 



In Fig. ini we plot a fc(x, /3) as a function of a; for various /3, and we confirm that our analysis agrees well 
with results from Brownian simulations. Comparing Fig. [3] with Fig. H] for f{x,y) shows that Tc{x,f3) is 
roughly by a factor (21na)^ faster compared to T{x,y). 



5 Discussion and conclusion 

We generalized here the narrow escape problem to a degenerated geometry defined by a flat cylindrical 
compartment of height h and radius R, where the absorbing hole is a small circular disk of radius a 
located centrally on the lower surface. We analyzed the problem for a laterally closed and open cylinder. 
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Figure 9: Conditional MFPT Tc{x, jS) to reach the hole in a laterally open cylinder for particles that are 
released on the surface opposed to the hole. The numerical values for fc(x,/3) are calculated using (|65|) 
and a = 50. The data points are the results of Brownian simulations with 10000 particles and a = 50. 



Because a uniform analytic expansion of the solution in the whole domain is not possible, we derived two 
different expansions in the two subregions il; and ilo, and then matched them at the boundary between 
the two subcompartments. 

We first analyzed the narrow escape time (NET) r(r, z) to reach the small hole in a laterally closed 
cylinder. For a flat cylinder with h <^ R and a, we obtained that the NET ([38|) is given by 



where the function ao{/3) E [0.07, 0.25] (Fig. [5^). Although was derived in the condition that h ^ R, 
our numerical results suggest that it remains valid until h ^ R (section I2.3.ip . In particular, for a 
cylinder with /i/a ^ In (■^) we recover the well known NET approximation r [36l [TOl [M], which 

was derived for non degenerated geometries (isoperimetric ratio of order 1 and no bottle neck). However, 
our analysis for the cylinder revealed that r = is already a very good approximation for an oblate 
volume with an isoperimetric ratio that can be very different from 1 (although h/a 3> In (■^), the ratio 
h/R can still be small). Formula can also be used to estimate the rate constant of a key chemical 
reaction during the early stage of phototransduction, which is the rate for diffusing cGMP molecules to 
reach the phosphodiesterase enzyme located on the surface of a narrow cylinder, located in the outer 
segment of a rod photoreceptor [23l [21] . 

In a next step, we used our method to analyze the narrow escape problem for a flat and laterally open 
cylinder, which is relevant for synaptic transmission. In many cases, the geometry of the synaptic cleft is 
well approximated by a laterally open cylinder [1] I22i 134) , where neurotransmitters are released into the 
synaptic cleft from the presynaptic terminal, located on the upper surface. The neurotransmitters move 
in the synaptic cleft (Fig. [1]) by Brownian diffusion, and they either activate receptors clustered in the 
postsynaptic density located on the postsynaptic terminal (corresponding to the lower cylinder surface) , 
or they leave the synaptic cleft through the lateral boundary without binding to a receptor. 

We estimated the conditional probability p{r, z) that a particle starting at position (r, z) reaches the 
small hole before leaving the cylinder. By identifying the small hole with the postsynaptic density where 
receptors are clustered, we estimated the fraction of released neurotransmitters that reach the receptor 
area before leaving the synaptic cleft. Using our analytic solution, we studied the impact of the synaptic 
cleft geometry as well as the location of neurotransmitter release. In Fig. [3 we plotted the probability 
to reach the postsynaptic density before leaving as a function of the release position for various cylinder 
height, and we found that it is very sensitive to the release position and the width of the cylinder. We 
conclude that, in order to achieve an efficient activation of postsynaptic receptors, the presynaptic and 
postsynaptic densities should be properly aligned such that the neurotransmitters are released opposite 




(68) 
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to the receptors. Finally, we computed the conditional mean time to reach the small hole before leaving 
through the lateral boundary (see (j65p V For a wide cylinder, we found that the conditional mean time 
is roughly by a factor (21nQ;)^ faster compared to the NET in a closed cylinder (see l\67\i ]. 

We shall now present some numerical estimates for neurotransmitters that need to activate receptors 
clustered in the postsynaptic density with a radius a — 50nm, when the synaptic cleft has a height 
h = 20nm and a total radius of i? = 500nm, so that a — R/a — 10 and (3 — h/a = 0.4. When the 
transmitter are released at distance r away from the center, the conditional probability is approximately 
given by truncating ()43|) at n = 



which is a very good approximation, as shown in Fig. [S] Using a diffusion constant D = 200/im^/s, we 
obtain for the mean times t(0, h) « 17/is and Tc(0, h) = \[is. 

The exact solution for the mean time and the conditional probability were obtained here by using 
the patching eigenfunction expansion approach, and from these approximations, we derived in the limit 
of large aspect ratio Rjh 3> 1 the asymptotic behavior. Our approach works well because of the radial 
symmetry due to the absorbing trap that is is located at the center of the cylinder. This situation accounts 
well for the postsynaptic density located at the center of the post-synaptic terminal. However, using our 
method it would be difhcult to treat the case of multiple non-concentric traps, and in this case a different 
approach based on matching asymptotic analysis should be more appropriate |21[ [6] . The analysis should 
start with an explicit representation of the Green's function for a cylinder. Once an inner solution is 
determined near each trap, it should be matched to the outer solution [3T]. This method should allow to 
study the effect of the trap positions and trap clustering on the synaptic current, which was only partially 
discussed in [331 HI] ■ 
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APPENDIX 

A Equation for the parameters in the Umit /3 <C 1 

In order to find the asymptotic equations for the coefficients a„ for a ^ 1 and /? — > 0, we introduce the 
scaled quantities 

Pn = PPn = triTTTT ' = ^ fcn s , 70 = P70 = 



In the limit /3 — > we have /„ — >■ cx) and fc„ — > oo for every n > 0, and the asymptotic behaviour of /3„ 
and ci„ is 



Using (j36|) . the asymptotic equation for the coefficients a„ in the limit /3 — > is 



oo ^ 

-CnO ■ (70) 

V 27r 



m=0 

Truncating (j70p at various levels n gives the approximations 

= ^rwO.lO, n = 

an 5 1 

-7= = 0.085, n=l 

V2 6 7r2 

^ ^ lIi-« 0.078, n = 2. 

V2 60 7r2 



Numerically we find from (j70p 



^ « 0.071 (71) 
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